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NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
TECHNICAL NOTE D-291 

LAMINAR BOUNDARY LAYER BEHIND A STRONG SHOCK MOVING INTO AIR 

By Harold. Mirels 


SUMMARY 

The laminar wall boundary layer behind a strong shock advancing into 
stationary air has been determined. Numerical results have been obtained 
for shock Mach numbers up to 14 using real gas values for density and 
viscosity and assuming Prandtl and lewis numbers of 0.72 and 1, respec- 
tively. The numerical results for shear and heat transfer agree, within 
4 percent, with a previously presented approximate analytical expression 
for these quantities. A slight modification of this expression results 
in agreement with the numerical data to within 2.5 percent. Analytical 
expressions for boundary-layer thickness and displacement thickness, cor- 
rect to within 4 percent for the present data, have also been obtained. 


INTRODUCTION 

Until 1957, all studies of the compressible laminar wall boundary 
layer behind a moving shock had assumed that the product of viscosity 
times density pp was constant throughout the boundary layer (e.g., 
refs. 1 to 5). This assumption becomes less valid as shock strength 
increases. Consequently a study was undertaken, at that time, to 
evaluate the effect of variable pp on the laminar boundary behind a 
strong shock moving into air. Real gas properties (ref. 6) were used. 

The boundary layer was assumed to be in thermodynamic equilibrium, and 
the Lewis number was assumed to be 1. Shock Mach numbers from 4 to 14 
were considered. 

The initial phase of the investigation was for a Prandtl number of 
1. Numerical results for shear and heat transfer were reported in ref- 
erence 7. Analytical expressions were used therein which correlated the 
shear and heat-transfer data to within 3 percent. The analytical expres- 
sions were generalized (by comparison with constant pp solutions) to 
account for Prandtl number not equal to 1, but the validity of this 
generalization was not established. 
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The second phase of the investigation was for a Prandtl number of 
0.72. Numerical results were obtained prior to the spring of 1958. 

These results were not published at that time in the hope that more ex- 
tensive data would be obtained. However, the press of other research 
projects prevented the continuation of the study. Recent informal dis- 
cussions with Dr. R. Hartunian, of the Cornell Aeronautical Laboratory, 
and Dr. N. Kemp, of the Avco Everett Research Laboratory, have indicated 
that the unpublished data are of current interest. In particular, the 
numerical data are of interest for correlating experimental wall heat- 
transfer measurements in shock tubes and for estimating the flow non- 
uniformity in low-density shock tubes. Hence the present publication of 
the data was undertaken. 

The present report includes the Prandtl number 1 data previously 
reported in reference 7, in order to make that data more generally avail- 
able. It is also shown herein that the approximate analytical expression 
for shear and heat transfer, developed in reference 7, is correct to 
within 4 percent for Prandtl number 0.72. Finally, approximate expres- 
sions are developed herein for boundary-layer thickness and displacement 
thickness which agree with the numerical results for Prandtl number 0. 72 
to within 4 percent. The latter expressions are useful for estimating 
flow nonuniformity in shock tubes (e.g., refs. 8 and 9). 

Numerical solutions for the wall laminar boundary layer behind strong 
shocks in oxygen recently have been presented in reference 10. Real gas 
properties, including the effects of Lewis number other than 1, are con- 
sidered. The correlation of theoretically derived heat transfer with ex- 
perimentally observed heat transfer in a shock tube (for both laminar and 
turbulent cases) is also discussed therein and in reference 11. In 
addition, reference 10 presents approximate analytical expressions for 
sheax and heat transfer which are similar to those developed in reference 
7 but which include the effect of Lewis number. 

The author is indebted to Richard J. Wisniewski for aid in obtaining 
the boundary-layer thickness correlations of the present report. 


ANALYSIS 

Consider the laminar boundary layer behind a shock moving into a 
stationary fluid. The problem is a steady-state one in a coordinate 
system fixed with respect to the shock (fig. l(a)). In this coordinate 
system the wall moves with the shock velocity = M s a i* 
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Equations of Motion 


If the boundary layer is assumed to be in thermodynamic equilibrium 
with a Lewis number of 1, the equations of motion for x > 0 are (e.g., 
ref. 12) 


dpQ + S( py) = o 

ox dy 



where H is the local stagnation enthalpy 


H = h + 



(la) 

( 11 ) 

(lc) 

(l d) 


( 2 ) 


The other symbols are defined in the appendix. From equation (la), a 
stream function exists such that 



Introducing new variables 



J 


( 3 ) 


(U 
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such that 


S d 5r] S 

ox c)£ Sx OT) 

s _ _p^_ j^t _s_ 

5y p w lf2xv w Stj 


and assuming ^ to have the form 



t = V2u e |v w f (r|) 

(6) 

reduces equations (lb) 

and (lc) to the form 



[Cf"] ' + ff" = 0 

(7a) 

B*] 

+ fg ' + IbH 1 - ■ 0 

(7b) 


with boundary conditions 

f(0) = 0 g(0) = ^ ^ u 2 

f ' (o) = 1 

f'(co) = : 

where 

C s 
U = 
g = 

In order to integrate equations (7), it is necessary to express G and 
a as functions of f and g. This is generally done by determining C 
and cr as functions of h/h w (for each free-stream pressure) and noting 
that 


g(oo) = 1 


( 8 ) 




V 


1. 

u^ 


( 9 ) 


H e H e + 2H e ^ > 


t 1 


cy p- 
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where 


fc w 


- (“ i?) w = ^ V5; 


y T 

H e 


p w J 2 t v v 

f u e 


[t] “ l( T] ) ] 


= 5^ = ^w 

u e 26 “p e 


■'w 


|26u. 



( 10 ) 

are 


f"(0) 

(11a) 

h'(0) 

(lib) 


(11c) 

- I(»)l 

(lid) 


121 


I(ti) 


■/’(>- a*. 


( 12 ) 


The integral l(t]) can be evaluated, after integrating equations ( 7 ), if 
p/p is known as a function of h/h . 


Integration of the Energy Equation for 0 = 1 
For 0=1, equations ( 7 ) and (8) become 


[cf] 

+ ff" = 0 

(l3a) 

[Cg’]' 

+ fg’ = 0 

(13b) 

f(0) = 0 

0? 

> 

II 

0 

'"5) 


f ' ( 0 ) = u 

1 — 1 

II 

IT 

*"So 

(13c) 

f*(») = 1 

j 
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Equation (l3b) can be integrated to yield 

g - g(o) = [frir ] (1 - « (0)) (i*> 

which is a form of the Grocco relation between velocity and stagnation 
enthalpy. If g(0) = 1, then 


g = 1 (15) 

so that the stagnation enthalpy is constant across the boundary layer. 
Recall that in a shock fixed coordinate system, stagnation enthalpy is 
conserved across the shock wave. It can then be shown that g(0) = 1 
provided h w = h^. The latter relation is approximately valid in the 

shock tube because the relatively high heat capacity and conductivity of 
the wall tends to maintain the wadi at its original temperature (e.g., 
ref. 7), and the pressure effect on enthalpy is small. When equation 
(15) applies, h’(0) / 0 and there is heat transfer to the wall (in con- 
tra distinction to the semi -infinite flat -plate problem where there is 
no heat transfer for the g = 1 case). 


Numerical Solutions 

Numerical solutions of equations (7) and (8) have been obtained 
which correspond to shocks of strength up to M s = 14 propagating into 

air at T]_ - 522° R and pj_ = 0.001 and 0.01 atmosphere. The Runge- 
Kutta method of integration was used with a step size At) = 0.02.^ The 
wall temperature was assumed to remain at 522° R for these calculations. 
Conditions behind the shock were determined by using thermodynamic charts 
for equilibrium' air (ref. 6). 

In the integration of equations (7) it was assumed that a was 
constant across the boundary layer. This is reasonable in view of the 


^This step size was established for the relatively low shock Mach 
number cases by decreasing the step size until the shear and heat trans- 
fer results were insensitive to a further change. Since decreases 

with increase in shock Mach number, a further check of the validity of 
this step size should have been made at M s = 14 - but this was not done. 
However, Dr. Richard Hartunian, in a private communication, has stated 
that he checked the effect of step size at M s = 14 and found that 
At) = 0.02 leads to accurate results for shear and heat transfer. He 
also stated that the computations of reference 10 were also made with 
At) = 0.02 (by error, ref. 10 notes the step size as being At) = 0.01). 
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rough computations of reference 13. The following analytic expression 
for C was used: 


= 1- 5 h ' 1 - °:1 481 + 0. 0028(t- - 

V h / h v h / h w ' h w 


5. 7 4X10“ 5 



(16) 


Equation (16) agrees to within 3 percent (for T w = T^ = 522° R and the 
specified ranges of pp and Mg) with computations hased on the thermo- 
dynamic charts for equilibrium air and the Sutherland viscosity- 
temperature relation. The thermodynamic charts provided p/p w ancl T/T w 

as a function of h/h w (for each specified T w and free-stream pressure) 
and C could then be found by using the Sutherland relation 


P 

Vw 


-j , 198.6 

/ t W/2 

T_ + 198. 6 VT J 

T w T w 


(17) 


The constants in equation (16) were chosen so that C has specified 
values at h/h w = 1, 10, 40 and has the correct slope at h/h w = 1. As 

a result, equation (16) provides an accurate representation of C for 
1 < h/h w < 42, which is the range required for the present investigation. 

Equation (16) tends to underestimate C for h/h w greater than 42. 
Similarly, p/p w can be approximated by 

0.02586 A + 0.94828 

-£- = A (18) 

P w A _ 0.02586 

h w 

for the range of pq and M s under consideration. The constants in 
equation (18) were chosen so that p/p w has specified values at 
h/h w = 1, 10 and has the correct slope at h/h w = 1. Equation (18) is 
correct to within a few percent for 1 < h/h w < 30 but is correct only 
to within 10 percent for h/h w near 42. It should be remembered that 
equation (18) is used only to evaluate l(q)- 

In order to determine the effects of Prandtl number, equations (7) 
were integrated for a = 1.0 and a = 0.72. To determine the effect 
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of variable C, equations ( 7 ) were integrated for C = 1 and for C as 
defined by equation (16). The results are summarized in table I. The 
results for C = 1 represent an extension of the results of references 
2 and 3 to values of U beyond 6. The shear and heat-transfer results 
for a = 1 and C = equation (16) have been reported in reference 7. 
The results for a = 0. 72 and G = equation (16) are new and represent 
the primary numerical contribution of the present paper. 


CORRELATION OF DATA 


Reference 7 proposed that the laminar -boundary-layer shear and 
heat transfer for variable C and constant a could be estimated 
from 

t”(P - ) = -0.489^1 + 1-665 U C?' 29 (19a) 

U — 1 ^ 


h' (0) 
h r - h w 


0.489-y/l + 1.665 U c °‘ 29 a ( 0 ’ 4:8+0 ' 022 U ) 


(19b) 


where h r is found from 


^ = 1 + (U - ir 


(0. 39-0. 023 U) 


2h £ 


(19c) 


A comparison of equations (l9a) and (l9b) with the numerical results ob- 
tained by integrating equations (7) is also given in table I. For 
o = 1, equations (19a) and (l9b) are correct to within 1 percent for C 
a constant and to within 3 percent for C varying according to equation 
(16). (The accuracy of eqs. (l9a) and (l9b), for o » 1, was previously 
pointed out in ref. 7. ) The present numerical solutions for o = 0. 72 
indicate that equations (19a) and (l9b) are correct to within about 4 per- 
cent for the data considered herein. 

Equations (19a) and (19b) would agree with the o - 0.72 data to 
within 2.5 percent if the factor were replaced by qO- 265 in 

these equations. The exponent 0.29 was proposed in reference 7 on the 
basis of the numerical solutions for the a = 1 cases. The use of the 
exponent 0.265 improves the correlation with the a = 0. 72 data, but 
at the expense of a poorer correlation of the a = 1 data. Since the 
a = 0. 72 solution is the more realistic one, for air, the exponent 
0.265 should be used in equation (19b) when making estimates of the 
wall heat transfer due to strong shocks moving through air. 

Reference 10 obtained numerical solutions for the shear and heat 
transfer associated with strong shocks moving through oxygen, assuming 
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a = 0.72 and different Lewis numbers. The data for Lewis number equal 
to 1 were correlated to within 1 percent by equations (19a) and (19b) 
with C e raised to the 0.24 power. The difference between the ex- 
ponents 0.265 and 0.24 is relatively small, considering that the 
viscosity law, the shock solutions, and the real gas properties used in 
reference 10 were different from those used herein. 

The displacement thickness &* in the steady coordinate system (fig. 
1(a)) is of interest for computing attenuation and nonuniformities m shock 
tubes (e.g., refs. 8 and 9). The displacement thickness is found from 
l(oo) and lim (f - T]) by substituting into equation (lid). Values for 

these quantities, obtained from the numerical integration of equations 
(7), are given in table I. The following approximate formula can be used 

to estimate 5*: 


where 


V u e Pe 
2xv w p w 


Cg' 37 lim (t ] - f) 
T] co 



3 


lim (,, - f)l s I - I54 < 1 - u ) 

r) -*oo _Ref 3 "Vl + 1- 022 U 


( 20 ) 


( 21 ) 


[l(°°) ] 


Ref. 3 


h-p h„ o ul 

— - — (U - l) 2 — 

1.154 h e h e 1.569 2h e 

-y/l + 1.022 U a 0. 47+0. 029 U -^/l + 0.993 U q0.045(U-l) 

( 22 ) 


Equations (2l) and (22) were obtained in reference 3 and correlate the 
constant property solutions to within 1 percent. The factor C^’ 57, in 

equation (20), corrects these expressions for the effect of variable C. 
The ratio of the value of &*, as obtained from equation (20), to the 
value of 5* obtained by numerical integration is given in table I. It 
is seen that, for the cases considered herein, equation (20) is correct to 
within 3 percent. 

The boundary-layer thickness 6 in the unsteady (wall stationary) 
coordinate system (fig. 1(b)) may be defined as the value of y at which 
the velocity relative to the wall reaches 99 percent of its free -stream 
value. That is, it corresponds to 
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or 



u - r 
u - 1 


0.99 




> 


f' = 0.99 + 0.01 U 


(23) 


This boundary-layer parameter is of interest, in experimental shock-tube 
studies, as a measure of the extent to which the boundary layer extends 
into the free stream. 


The value of rj corresponding to y = 6 is denoted herein by Tjg. 

The boundary-layer thickness 5 can then be computed from equation (lie) 
by using the values of rjg and l(oo) given in table I. (More properly, 

l(ir]g) should be used when computing 6 from equation (lie), but the use 

of l(co) should introduce a negligible error. ) The resulting values for 
5 are included in table I. 


The following approximate formula can be used to estimate 5: 


V u e 
2xv w P w 


<t 48 [n 6 - K-5] 


Ref. 3 


(24) 


where 


&4tof. 3 


5.20 

VI + 0.543 U 


(25) 


and l(co) is found from equation (22). (Again, l(r]) has been evaluated, 
for convenience, at r\ = <» rather than T]g. The error should be unim- 
portant. ) Equation (25) was found by assuming rjg to be of the form 

A/^/l + BU and determining A and B from the U = 1 and U = 6 
constant property results of reference 3. (Eq. (25) differs slightly from 
a similar formula given in appendix D of ref. Q and is the more accurate 
expression if the edge of the boundary layer is taken to be 

u/u e = 0.990.) The coefficient C^ 48 , in equation (24), is the approxi- 
mate correction to account for variable C. The ratio of 5, as obtained 
from equation (24), to the value obtained from the numerical integration 
of equations (7) is given in table I. Equation (24) is correct to 4 per- 
cent for the cases noted herein. 


tr. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, December 7, 1960 
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APPENDIX - SYMBOLS 
a local speed, of sound 

c pp/(pp) w 

f function of T) defined By eq. (6) 

g H/H e 

H stagnation enthalpy in shock stationary coordinate system 

h static enthalpy 

1(0 integral defined by eq. (12) 

M s Mach number of the shock 

p pressure 

q w heat-transfer rate 

R gas constant 

T absolute static temperature 

u Uw/ U e 

u,v velocities parallel to x, y axis 

u,v velocities parallel to x,y axis 

u s velocity of shock in x, y coordinate system 
u^ velocity of wall in x,y coordinate system 

x, y steady coordinate system (fig. 1(a)) 

x, y unsteady coordinate system (fig. l(b)) 

z molecular weight ratio 

B fluid velocity - boundary-layer thickness (value of y corres- 

ponding to u/u e = (u - u w )/(u e - u w ) = 0.99) 

B* fluid-boundary-layer displacement thickness in steady coordinate 



system, 
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rj similarity parameter, eq. (4) 

^5 value of r\ corresponding to y = 6 

M- coefficient of viscosity 

v kinematic viscosity 

£ function defined by eq. (4) 

p mass density 

0 Prandtl number 

T w local shear stress exerted by fluid on vail 
ijr stream function, eq. (3) 

Subscripts: 

e flow external to fluid boundary layer 

r quantity evaluated for zero heat transfer 

w quantity evaluated at the vail 

1 undisturbed flow ahead of shock 

Superscript : 

' denotes differentiation with respect to r\ 
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(b) Unsteady coordinate system. 

Figure 1. - Coordinate systems used to study boundary layer behind a 
shock wave advancing into a stationary fluid. 
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